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Abstract 



p\ ' We show that the 2-inatrix string model corresponds to a coupled system of 2 + 1- 

j^ ■ dimensional KP and modified KP ((m)KP2+i) integrable equations subject to a spe- 

cific "symmetry" constraint. The latter together with the Miura-Konopelchenko map 
for (m)KP2+i are the continuum incarnation of the matrix string equation. The 
(in)KP2+i Miura and Backlund transformations are natural consequences of the un- 
derlying lattice structure. The constrained (m)KP2+i system is equivalent to a 1 -|- 1- 
dimensional generalized KP-KdV hierarchy related to graded SL(3, 1). We provide an 
explicit representation of this hierarchy, including the associated W(2, l)-algebra of 
the second Hamiltonian structure, in terms of free currents. 
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1. Introduction 

One of the remarkable features of integrable Toda lattice hierarchies, describing the multi- 
matrix string models [|l|, is the existence of a natural elegant way p[ to extract from them 
continuum differential integrable systems without taking any continuum (double-scaling) 
limit. These systems are of generalized KP-KdV type and are associated with one fixed lattice 
site. This proves to be an efficient approach for exact calculations of multi-matrix model 
correlation functions and exhibits deep relations with topological field theory. In particular, 
the two-matrix model in the formalism of provides a systematic unified framework for 
exact treatment via integrable hierarchies of c = 1 string theory (for related approaches, cf. 
0]). The principal instrument are the Wi+oo constraints on the two- matrix model partition 
function, which are the facets of the pertinent string equation. 

The aim of the present note is to construct a differential integrable system out of the 
Toda lattice hierarchy, underlying the two-matrix model, which directly incorporates in itself 
the whole information of the string equation, rather than imposing it as an infinite number 
of constraints on the discrete model partition function. We achieve this by assuming that at 
least one of the matrix potentials in the two- matrix model is of finite order p. In the simplest 
nontrivial p = 3 case the continuum system, we are looking for, turns out to be a coupled 
system of 2 + 1-dimensional KP and modified KP ((m)KP2+i) integrable equations subject 
to two additional constraints. The latter are identified as the known Miura-Konopelchenko 
map 0] and a specific KP "symmetry" constraint |^. We demonstrate that these two 
differential constraints embody the matrix string equation. In particular, the two-matrix 
model "succeptibility" satisfies the above string-constrained KP2+1 equation. 

Lattice translations between different sites in the context of Toda hierarchy are shown 
to generate the well-known (m)KP2+i Backlund transformations @, |^ which are canoni- 
cal symmetry of the (m)KP2+i Hamiltonian system. Furthermore, from the constrained 
2 -|- 1-dimensional (m)KP system we derive an equivalent unconstrained 1 -|- 1-dimensional 
integrable system of generalized KP-KdV type related to graded SL(3, 1) algebra 0, 0. 
Systems of this type were previously obtained in ||10|, [Tl| by imposing ad hoc second-class 
constraints on multi-boson KP hierarchies. The corresponding second Hamiltonian structure 



is the nonlinear W(2, 1) algebra, generalizing the Wq^ algebra [12|. We provide an explicit 



free-field representation of W(2, 1) (cf. also [[13 , 



2. Discrete Integrable Systems from Two-Matrix String Models 

This section starts by collecting some basic results on the discrete integrable systems - 
generalized Toda-like lattice systems, which are associated with multi-matrix string models 
0. The associated Toda matrices contain in general finite number of non-zero diagonals 
unlike the case of the full generalized Toda lattice hierarchy [|T^]. We present the relevant 
solutions for the Toda matrices in a form appropriate for establishing connection with 2 + 1- 
dimensional integrable models. 

Specifically, we shall consider the two-matrix model with partition function : 

Zj^[t, i, g]= dMidM2 exp - <^ ^ i, Tr M[ + ^ t, Tr M^ + gTi M1M2 \ (1) 

•^ It-=1 s=1 J 



where Mi 2 are Hermitian Nx N matrices, and the orders of the matrix "potentials" pi^2 niay 
be finite or infinite. In ref.p] it was shown that, by using the method of generahzed orthog- 
onal polynomials |T^, the partition function (|1|) and its derivatives w.r.t. the parameters 
(tr, is, g] can be explicitly expressed in terms of solutions to generalized Toda-hke lattice sys- 
tems associated with (Q). The corresponding linear problem and Lax (or "zero- cur vat ure" ) 
representation of the latter read M : 
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The subscripts — /+ denote lower/upper triangular parts, whereas (+)/(—) denote up- 
per/lower triangular plus diagonal parts. The parametrization of the matrices Q and Q 
is as follows: 
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From now on we shall consider the first evolution parameters ti,ti as coordinates of a 
two-dimensional space, i.e., ti = x and ti = y, so all modes ak{n),bk{n) and Rn depend on 
{x,y;t2, . . .,tp^;t2,. ■ ■ ^tp^j ■ 

The second lattice equation of motion (^) for s = 1, using parametrization (j^), gives : 



^x^n 



d^ 



hk{n) 



Rn (boin) - boin - 1)) 
bk^iip) ~bk+iin- 1) 



dMn) = b,{n) ~ biin - 1) (9) 
k>2 (10) 



\Rn+l ■ ■ ■ Rn+k J Rn+1 ■ ■ ■ Rn+k 

Similarly, the first lattice equation of motion (|^) for r = 1 gives : 

dybo{n) = Rn+i - Rn , dybk{n) = Rn+ibk^i{n + 1) - Rn+kbk-i{n) 



(11) 



for k > 1. From the above equations one can express all bk{n ± i) , k > 2 and Rn±e (^ ^ 
arbitrary integer) as functionals of bo{n), bi{n) at a fixed lattice site n. 

In complete analogy, the lattice equations of motion (H) for r = 1, s = 1 read explicitly: 



dxa-oin) = Rn+i — Rn 



dxttkin) = Rn-k+io-k-iin) - Rnak-i{n - 1) 



(12) 



(with A; > 1) and 

dyao{n) = ai(?2 + 1) — ai{n) 
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with A; > 1. Following (|), (0), (||) we obtain the "duahty" relations: 

dyhi{n) = dxRn+i , d:^ao{n) = dybo{n) , d^aiin) = dyRn 
From the above one gets the two-dimensional Toda lattice equation: 

dy In Rn = ao{n) - ao{n - 1) -* 9^5^ In Rn = Rn+i - 2-R„ + Rn-i 



(13) 



(14) 



(15) 



Eqs.(|T2|)-(|l^) allow to express all ak{n ± €) , k > 1 and Rn±e as functionals of ao(^) 
and Rn (or ai{n) instead) at a fixed lattice site n. Furthermore, due to eqs. (|T^) and ([TT|) for 
fc = 1, all matrix elements of Q and Q are functionals of bQ{n), bi{n) at a fixed lattice site n. 
Alternatively, due to (|14D we can consider ao{n) and Rn+i as independent functions instead 
oibo{n),bi{n). 

Let us also add the explicit expressions for the fiow eqs. for bo{n) , bi{n) , Rn+i resulting 
from M) and (IBD : 
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3. Solution of String Equation for Finite-Order Matrix Potentials 

From now on we assume that one of the matrix potentials in (|l|), e.g., the second one has 
a finite order p2, which implies that the matrix Q has a finite number of diagonals below 
the main diagonal (cf. (|^)). In this case the lattice equations of motion impose additional 
constraints relating the two independent functions 6o(^) and bi{n) (or -R„+i). 

More precisely, we find that for p2 < oo the lower triangular plus diagonal part of Q is 
expressed through Q on the space of solutions to (p!^)-(p!5|) in the following form (in what 
follows, we shall ignore arbitrary integration constants without loss of generality) : 

Q = gf!T' + (-x) 1 



'(-) 



.9 



where g is the two-matrix model coupling parameter appearing in the string equation 
and I+nm — ^n+i,m ■ In particular, for the lowest non-zero diagonal of Q one has: 



-Rn 



ap2-i[n) - Kn- ■■ Rn-p2+2 - Q'^n,n-(p2-l) (^9) 

which follows from lattice eq. (|13D for k = p2 — ^ and from the explicit parametrization of Q 



Eq.(|T^) can be proved by induction w.r.t. k = p2 — I, ■ ■ . ,1,0 starting from (0), upon 
comparing eqs.(|l2D-(|l3]) with the component form of the matrix equations of motion dyQ"^ = 
Q{+) , Q^ and dxQ'^ = Q* , Q- for arbitrary integer power s. 

The remaining lattice equations of motion - the first equations of (|l^) and (|15|), imply 
the following two additional constraints on the independent functions bQ{n) and Rn+i (or 
hin)) 

dyboin) = - + d, (QZ-') , dyRn+i = d, {Q':,Xi,n) (20) 

In fact, eqs. (pO]) are nothing but the component form of the string eq.(^. Indeed, it can 
be rewritten in the form (upon using ([T8|) and (|D,(H)) : (dy — d/dip^A Q = -1. , i.e. 

d Id d 

-bo{n) = dybo{n) , -j R„+i = dyRn+i , —f bk{n) = dybk{n) (21) 



for k > 1. Now, inserting ( |T6|) into (pO|) we find that the latter two equations precisely 
coincide with the nn and n+1, n component of the string eq. ([21|) . In particular, the evolution 
parameter tp2_i is not independent but rather essentially coincides with y = ti. It turns out 
in what follows that it is more convenient to use y = tp^-i as a second space coordinate 
instead of y. 

To conclude this section, let us note that there is a complete duality under pi < — *• 
P2 when the order pi of the first matrix potential in (|l|) is also finite. We obtain the 
exact analogues of eqs.(|l8D, (|20|) and (plf) by interchanging pi < — > p2 , x = ti < — > ti = 

y , ip2-i ' — ' ^Pi-i > Q{~) ' — ' Q{+)- 

4. Derivation of mKP2+i and KP2+1 from Lattice Systems 

Using the parametrization (0)^(11), the equations of the auxiliary linear Lax problem (0),(^ 
acquire the form: 

P2-1 

A^n = i^n+i + ao(n)?/'„ + Y^ ak{n)i)n-k (22) 

fc=i 

1 d ^^~^ bk(n) 

aV^" = Rni'n-l + bo{n)i)n + "Y — V'n+fc (23) 

9 (^A j^^-^^ -K„+l . . . Kn+k 

dxipn = -Rn'ipn-l , dylpn = 1pn+l + ao{n)^n (24) 

In particular, one has from the first eq. (|2^) : 



(-1)'= 1 1 

ipn-k = s d^^ d^ . . . -—d^ipn (25) 

-Kn-fc+l ^n-k+2 ^Mi 

Inserting (|2^-(P3D into ( P^ one obtains the following two-dimensional differential Lax spec- 
tral equation for ip^ at a fixed lattice site n ^ : 

Lp,-i{n)i Jn = (26) 

^In eq.(|26|) we have hidden the spectral parameter A by using the identity -^r^ — ipn = dytpn = {dy — A) t/jn 
, which foUows from (|1),(|) and ® . 
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where in the second hne of eq.(|27|) we have used the lattice equations of motion (0), i.e., 
Q- = [Q^^^^] , to express the coefficient functions of the 2-dimensional Lax differential 

operator in terms of the independent 60,1 ('^)- 

Similarly, the lattice flow equations w.r.t. ig (^) can be equivalently written as differential 
flow equations: 
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^l,^ = -M,{n)i,^ , M,{n)^{-iyd: + Y.{-lffi\ho{n)M{n))d': (28) 

with s = 2, . . . ,p2- Explicitly we have for the coefficient functions: 

/a = sb,{n) , fi% = shin) + ('] [blin] - dMn)) , etc. 



(29) 



In particular, comparing (|28|)-(p9D and (p6D -(^7D, we note that the Lax spectral equation 
(p6|) is in fact the {p2 — l)-th flow equation (|28|). 
The compatibility conditions for (^) : 

d d y 1 

-^Ms{n) - ^^M.,{n) + M,.{n), Ms{n) =0 , r,s = 2,...p2 (30) 



dt 



dt. 



yield for r = p2 ~ 1 (i-e., tp2_i = y) a. system of 2 + 1-dimensional integrable nonlinear 
evolution equations for &o(^)) biiji) subject to the additional string-equation constraints (pO|). 
Let us consider in more detail the simplest nontrivial case P2 = 3 where: 

M2{n) = dl~2h^{n)dx . M:,{n) = -&i + ?,ho{n)dl- 2,[h{n) + hl{n) - dMn)) d. (31) 

Eq.(^0|) with r = 2, s = 3 yields firstly the non-dynamical equation: 

dM'^) = 9. (2&i(n) - dMn) + bl{n)) (32) 

which coincides with the first string eq.(pOD (for p2 = 3). Secondly, it yields an evolution 
equation for hQ{n) : 



4^feo(^) 

oh 
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dxao{n) = dMn) , ao{n) = ao{n) x = 2bi{n) - dxbo{n) + bl{n) (34) 

where we have used (|T^ and (P^). Combining (^2p with ( PB| ) yields an evolution equation 
for bi{n) : 



4ifei(n) = 9^61 (n) + 69,6? (n) + 3dyRn+i 
oh 



dxRn+i = dMn) 



(35) 



Finally, the second string eq-PDj) (for p2 = 3) yields: 

dyRn+i = d^ [dr,Rn+i + 26o(ri)i?„+i] = M*(n)i?„+i (36) 



where M2{n) denotes the operator adjoint to M2{n) in (|3lD . 

Now, it is straightforward to recognize eq.(^) as the 2 + 1-diniensional KP equation 
(KP2+1), eq.(0) - as the modified KP equation (mKP2+i), and eq.(|32D - as the associated 
2 + l-dimensional Miura-Konopelchenko map @] relating them. However, the present two- 
matrix model derivation reveals an additional constraint ( p6D relating KP2+1 and mKP2+i 
equations. This constraint can be identified as a particular example of a "symmetry con- 
straint" P within the framework of the standard 1 + l-dimensional KP hierarchy (KPi+i) 
formalism. The KP reduction, triggered by (|36|), is studied in section 6. 

As demonstrated above, both the (m)KP2+i eqs. (PB| ) and (P^D, as well as the Miura- 
Konopelchenko map ( p2D and the "symmetry constraint" (|36D KP reduction naturally arise 
from the underlying Toda-like lattice integrable structure associated with the two-matrix 
string model. Let us emphasize that the Miura-Konopelchenko map (0) and the "symmetry 
constraint" (|36|) explicitly embody the "string equation" (^) of the two-matrix model in the 
case P2 = 3. Thus, they are the 2 -|- l-dimensional analog of the Wi+00-constraints on the 
partition function (|l|) in the formalism of refs.0. 

Through the method of orthogonal polynomials the two-matrix model partition function 
(H) is given by: 

Z M[t,i,g\= const N\Whn , K = exp bo{n) (37) 

n=0 "^ 

Thus, its calculation reduces to finding the appropriate solutions of (m)KP2+i equations 
(0),(^) subject to the "string equation" constraint (p6|). In particular, as follows from ( p7[) 
and d^), the "succeptibility" dlln Zj^ = bi{N — 1) satisfies the string-constrained KP2+1 
equation (|35|),(|36|). 

5. Syminetries of Constrained KP2+1 and niKP2+i 

In the previous section we represented the differential KP2+1 and mKP2+i integrable sys- 
tems explicitly in terms of objects from the underlying generalized Toda-like lattice structure. 
Now we are going to show that this underlying lattice structure naturally exhibits the discrete 
Backlund transformations for KP2+1 and mKP2+i. 

On the lattice itself we have two obvious discrete symmetry operations: 

(a) lattice site translation n — > n—1, in particular, for the discrete Lax "wave" function 

1pn+l ^ i'n = -Rnlldxi^n+1 (38) 

(b) similarity (phase) transformations on the same lattice site, e.g., gauge transformations 

'ipn ^ V'n = e'^-lpn = $„V^n (39) 

First, we consider type (a) lattice symmetry (^). Going back to mKP2+i (0) and 
KP2+1 (p5|), we observe that bo{n) and bi{n) do satisfy them for any fixed site n on the 



6 



underlying lattice, i.e., if bQ{n),bi{n) satisfy (P5|),(P^), so do bo{n — l),6i(n — 1) . On the 
other hand, from the lattice equations of motion (j^), ([TT|) we obtain the following simple 
relations for the lattice shifts of the above functions: 



bl-\n) 



boin 

R-n - 



■1) = 
Rn+1 
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dybo{n 
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) , b[-\n) = b,{n-l) 
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b^{n)-^Mn)m 
(41) 



which can be viewed as discrete symmetry (Backlund) transformations for (m)KP2+i ( pSJ) 
and (^5|) . Obviously, type (a) Backlund transformation (^Of) preserves the (m)KP2+i "sym- 
metry" constraint (|36|). Note that the two- matrix model coupling constant g emerges here 
as a free Backlund parameter. 

Let us also note, that the Backlund transformation for bi{n) (^)) (taking into account 



(32)) is a canonical transformation for the 2 + 1-dimensional KP Hamiltonian system (cf. 
T^), i.e., it leaves invariant the canonical Poisson brackets and the Hamiltonian (here there 



is no difference between y and y) : 

{b^in){x,y) , b^{n)ix',y')} = ^dj{x - x')6{y - y') = {b^{n + l){x,y) , b^{n + l)ix',y')}m 

lo 

H{n) =2Jdxdy [{dMn)f - 46? (n) - ml^^] = H{n + 1) , d^Rn+i = dybi{n) (43) 
In the framework of the standard KPi+i hierarchy formalism with KPi+i Lax operators: 
L^KP = d^ + v + Y.^kd^^ ' LKP = d^ + J2ukd-^ (44) 



fc>i 



fc>i 



eqs.(^) have already been obtained [^ |^ (identifying v = bo{n), Ui = 2u = 2bi{n) ) via the 
method of pseudo-differential operator gauge transformations and are known as the first type 
of auto-Backlund transformations for conventional (m)KP2+i. However, in the present two- 
matrix model realization of mKP2+i and KP2+1 there is an additional important property 
of this auto-Backlund transformation — the additional constraint (^) expressing $ = $„ as 
a specific explicit functional of the (m)KP2+i functions u, v . This additional constraint on 
the auto-Backlund transformation results from the (m)KP2+i "string-equation" constraint 
(^). Both constraints (^) and (^T]) are absent in the general (m)KP2+i setting. 

Next, we consider type (b) lattice symmetry (|39|). It must preserve the form of the 
2-dimensional Lax and the corresponding tg-fiow generating operators: 



$: 



'A 



Ms{n) ($, 
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dt 



^ + Ms{n) , s = 2,...p2 
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M,{n)^{-iyd: + Y.{-^tfi'^(hoM)dl , b 



oW 



60 (n) - d^ln^n , bi{n) 



(45) 
6i(n)(46) 



fe=i 



where the last two equations represent type (b) Backlund transformation for unconstrained 
(m)KP2+i P, 0]. Comparing (^5|) with the similarity relation between the operators Ms{n) 
and Ms{n + 1) implied by type (a) lattice symmetry (pSj) : 
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-dx 



■n+l 
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^ + Msin + 1) 



d-\Rn+r) = ^ + Ms{n) , s = 2, 
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we find that $„ in (|i3|)-(^BD is again given by eq.(|T). However, with such $„ type (b) 
Backlund transformation (|46|) does noi preserve the string-equation constraint (|36|). There- 
fore, only type (a) (lattice translation) Backlund transformation (^Oj) survives in the two- 
matrix-model realization of (m)KP2+i. 

Concluding this section, let us consider also the continuum symmetries of the con- 
strained (m)KP2+i. In this respect we note further important differences between the 
string matrix-model realization of mKP2+i and KP2+1 integrable systems, on one hand, 
and the standard KPi+i integrable hierarchy, on the other hand. In the present 2 + 1- 
dimensional matrix-model formalism the number p2 of the usual KP time-evolution param- 
eters (ti = x,t2,... ,tp2_i = y,tp^) is finite, e.g., p2 = 3 in the above analysis, unlike the 
KPi+i case. Furthermore, besides them there is also an additional set of time-evolution 
parameters ti = y,t2, ■ ■ ■ ,tp^, which might be infinite in number, i.e., pi — >■ 00. Their 
corresponding flows, given by ([T7|), are both commuting with the tg-fiows (given by (|16])), as 
well as commuting among themselves. So the latter correspond to an (infinite-dimensional) 
Abelian symmetry algebra for the constrained string-matrix-model's (m)KP2+i system. 

The nature of the continuum symmetries corresponding to the t^-Aows could be better 
understood if one considers the constrained (m)KP2+i in terms of the equivalent 1 + 1 
dimensional integrable system (eqs. (^8])-(^) below). 



6. Equivalent 1 + l-Dimensional Formulation of Constrained 

(m)KP2+i 



The string-constrained (ni)KP2+i system (^)-(|36|), which describes the two-matrix model 
(for P2 = 3), can be explicitly reduced to an equivalent 1 + 1-dimensional generalized KP- 
KdV integrable system. Upon excluding in (^) ao(^) via (|3^) and (|32D , and substituting 
(P5|) into (^), the corresponding reduced 1 + 1-dimensional integrable system can be written 
in the following form: 



ots 

4|-6i(n) 
ots 






(48) 



dlRn+i + 3 {hi{n) + hl{nfj R^+i + 3bo{n)dx:Rn^ 

dXin) + blin) + 3b,{n)boin) + |i?„+i - h, (61(72) + bUn))] (49) 



dx dlbi{n) + Qb\{n) + 3dxRn+i + 66o(^)-Rn+i 



(50) 



Similarly, the ^2 flow eqs. acquire the form: 
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dxRn+i + 2bo{n)Rn+i 
2bi{n) + bl{n) - dMn) 
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—^bi{n) = d^Rn+i 
dt2 



(51) 
(52) 



In complete analogy with the 2 + 1-dimensional type (a) Backlund transformation (]40D , 
from the lattice equations of motion one immediately obtains the Backlund transformation 



for the 1 + 1- dimensional system (^S|)-(^) resulting from negative lattice site translation: 

bi-\n) = boin - 1) = boin) - 9,ln(i?„+i ---8, (26i(n) + bUn) - dM^^))) (53) 

b[~\n) = b^in - 1) = bi{n) - dMn) (54) 

Rn = Rn+1 - - - 5x (26i(n) + bl{n) - dMn)) (55) 



K. 



(-) 

n+l 



Again, as in (^Of ), the two- matrix model coupling constant g appears as a free Backlund 
parameter. 

The system (|48| ) - ([50|) (or (|5lD -(p2|)) possesses Lax representation which is 1+1-dimensional 
analog of (©-(i^) (with p2 = 3) : 
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-d-\Rn+r) + QZ:' + M,,^, 



n] 



L2in) = -^-'{Rr.+r) + 2b^{n) + {^,-boin)y 



for p2 



(56) 

(57) 
(58) 



In particular, one reproduces the lattice flow eqs. ([T6|) as compatibility conditions for the 
linear auxiliary Lax problem (|5B|). The systems (^5[)-(p]|) and (^T])-(^) themselves arise as 
compatibility conditions: 



dts 



M^{n),L2{n) 







dt2 



M2(n),L2(n) 







(59) 



with L2{n) as in (|58D, and M2,3(n) as in (|31|). 

Let us concentrate on the Lax operator (^) (indication of the associated lattice site 
n is suppressed for brevity) : Z2 = —D^^ R + 2bi + b^ + d^bo — D (26o) + -D^. The Hamil- 
tonian structures inherent in the integrable system (^8l)-(^) (or (^Il)-(^)) can naturally 
be identified using the general i?-matrix scheme of Adler-Kostant-Symes/Reyman-Semenov- 
Tian-Shansky [ll7| , |l^ in the context of the algebra of all pseudo-differential operators ^T>0 
. Namely, the first and the second Hamiltonian structures are given in terms of the following 
i?-matrix Poisson brackets: 



{(L|X),(L|F)}« 
{(L|X),(L|F)}g) 



■{L\\[RX,Y] + \[X,RY]) 
-{L\\[R{XL),Y] + \[X,R{LY)]) 



(60) 
(61) 



where L, X, Y are arbitrary elements of \I/P(9: L = Y^k>-oo D^ ^fc ^^^d X = J2k>-oo ^kD^ 
and < "I" > denotes the standard bilinear pairing via the Adler trace. {L\X) = Tr^ (LX) = 
J dxKes (LX). Both Hamiltonian structures are compatible provided R is antisymmet- 
ric operator on '^T>0 w.r.t. the Adler trace. There are three different natural i?-matrix 

^To avoid confusion, from here on D will denote the differential operator w.r.t. x, whereas dxf will 
denote derivative acting on a function. Note also the inverse ordering of differential operators and coefficient 
functions in the definiton of the Lax operators. 



operators on ^VO labelled by £ = 0, 1, 2 |g : 

oo 

Re = P>i — P<£-i , P>eX = 2_^ X^D , P<£_iX = 2^ X^^D (62) 



fc>£ 



k=e-l 



Now, one observes (cf. ||2^, ^) that the space of Lax operators L2 {p2 = 3), and more 
generally - the Lax operators for arbitrary (finite) p2 (|5l 



P2- 



L,,^, = -D-'R+J2{-lfD''f, 



k-nk p{P2-l) 



'_X)P2-l£)P5 



-1 



(63) 



fc=0 



with f^^~ simply expressed in terms of f^^~ from (|29| ) for s = p2 "~ 1? span a closed Ri- 
coadjoint orbit on ^VO with Ri as in (^). The corresponding first Hamiltonian structure is 
given by the i?-Kirillov-Kostant bracket (|60|) with L = L2 and R = Ri = P>i — P<o. The ^3- 
and t2-flow eqs.(P5|)-(|5(]D and (pT|)-(|52D are Hamiltonian equations of motion w.r.t. the first 
bracket (|60|) , which reads exphcitly { R{x) , hQ^x') } = —dx5{x — x') and { &i(x) , 61 (x') } = 
— \dx5{x — x'), and with HamiltoniansQ H^ = —j^TtaL2 for A^ = |, 2 . 

Concerning the construction of the second Hamiltonian structure, we note that Ri from 
(|6^) is not antisymmetric. However, there exists a symplectic mapping |20, 21, 0] from the 
i?i-brackets to the standard -Ro-brackets of the first KP Hamiltonian structure where now 
Rq ( |62D is already antisymmetric: 



{ £= 



X 



u 



^'))>{ 



^(0) 



X 



L 



(0) 



^)) 



(1) 



^(0) 



I''«L2eI''° (64) 



with L2 = — (D + bo)^ R + 2bi + D^. Therefore, we can now use the general formula 
(|6l|) to construct the second Poisson brackets compatible with the first bracket, however, 
with a caution. Namely, let us emphasize the following important point. Eq.(|6lD is given 
for generic pseudo-differential Lax operators L. Restriction of L to arbitrary submanifold 
does not necessarily lead to a consistent restriction of the corresponding i?-matrix Poisson 
brackets (|6TD (see, e.g. [^). Thus, we have to prove that the restriction of generic L to 



^(0) 



ID is in fact a consistent Poisson reduction. This can be done by adapting the proof 
in [|13| for the multi-boson reductions of the ordinary KP hierarchy. The final result for the 
second Hamiltonian structure for Lg reads Q : 



f{0) 



X 



f(o) 



^■>} 



(2) 



Tr. 



Zfx)^^zrV 



+ - 



dxResl 



?(0) 

-^2 



X 



xlf 



9-^Res( 



>o 



Ylf 



?(o) 
-^2 



Y 



(65) 



and the specific expressions for the nonlinear Poisson brackets among the coefficient fields 



can be found in refs. iTOl, ITl 



^The overall minus sign in the definition of Hj^ is due to the inverse ordering of ZJ's and coefficient 
functions in L2 ■ 

^The last term on the r.h.s. of ( |65[ ) is a standard Dirac-bracket term due to the absence of next-to-leading 

(0 

2 



differential order term in L, ( |64| ) 
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It follows from the general considerations in [TT|, that the coefficient fields i?, Bq, bi gener- 



ate the W(2, l)-algebra, which is a generalization of the nonlinear Woo-algebra [0. On the 
other hand, in |T3| we succeeded to represent the usual multi-boson KP hierarchies in terms 
of canonical pairs of free fields abelianizing the second Hamiltonian structure of ordinary 
KP hierarchy. Using the latter result, we find a similar result for L2 - explicit expressions 
of the coefficient fields (Rn+i,bo{n),bi{n)) in terms of free fields (currents) (ci, 61,02) : 

R= {d^ + ei + ci + 2C2) (<9^. + ci) ei 
&i = - {d^ + ci) ei - - {dx + C2) C2 , 60 = ei + ci + C2 (66) 

{ci{x) , ei{x')} = -dj{x - x') , {c2ix) , C2{x')} = -dj{x-x') , rest = (67) 

In turn, formulas (0),(0) provide an explicit free-field representation of W(2, l)-algebra 
embodied in the second bracket structure (p3|). Let us note, that (|5BD-(|5^) can also be 
obtained from the free current's realizations of the modified full KP hierarchy P] by a 
specific Dirac constraint reduction. Similar free- field representations exist for W(p2 — 1,1) 
Poisson-bracket algebras for the higher Lax operators (|63|) which describe generalized graded 
SL(p2, 1) KP-KdV integrable systems [|, |]. 

We conclude with a remark about the infinite-dimensional Abelian symmetry algebra of 
the 1 + 1-dimensional integrable system (^)-(l5^), (|58D -(p9|), generated by the t^-fiows ([T7|) . 
Inserting the solution for Q ([TB|) into ([T7|), we find that the corresponding vector fields -S- 
are Hamiltonian ones, with Hamiltonian functions T-Cr w.r.t. the first Poisson bracket (|60| 



of the following form: 7i,, = Hr+i = — ;^TrA-^2^^ for integer r > 1. The Lax operator L2 
is of the same form as (|57|)-(p8|), but with ao{n) = Q^„ + -a; as a zero-order term instead of 
Ql^^2b,{n) + bl{n)-dMn) . 

Outlook It would be interesting to extend the analysis of sections 4-6 beyond the simplest 
case of P2 = 3 and for higher multi-matrix models. The general case of Lax operators I/p2-i 
(|53|) with arbitrary finite differential part of order p2 — 1 is studied in ref . ||^ . From eq. (|5^) 
we expect then to find explicitly the generalized SL{p2, g)-KdV hierarchies with q restricted 
to g = 1. One recalls at this point that the generalized Kontsevich models |^ provide 
description of the continuous models (based on c < 1 minimal conformal models) only for 
g = 1 in the {p, g)-series and it would therefore be natural to study a possible relation. 
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